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Abstract— One of the main challenges in through-the- 
wall radar imaging (TWRI) is the long data acquisition 
and/or data processing time. The movement of the target 
during the long data collection time will cause smearing 
and target displacement in the through-the-wall image. 
Rapid data processing is critical to the real-time 
monitoring of targets behind the wall. In this paper, we 
present three techniques for fast data acquisition and 
real-time processing for TWRI. The first technique is 
based on compressive sensing incorporating the layered 
media Green’s function and is capable of dealing with 
targets behind multilayered building walls. Next, we 
formulate the beamforming algorithm for through-the- 
wall imaging with multiple-input multiple-output 
(MIMO) radar which is real-time in data collection. 
Finally, we propose a real-time processing algorithm for 
radar imaging through single or multilayered walls. 
Numerical results are presented to show the efficiency 
and the promising applications of the above techniques 
for fast data acquisition and real-time processing in 
TWRI. 


Index Terms— TWRI, fast data acquisition, real-time 
processing, compressive’ sensing, MIMO radar, 
diffraction tomography. 


I. INTRODUCTION 


The high-resolution and noninvasive imaging of targets 
through visually opaque obstacles, such as walls and doors 
has sparked a growing interest in through-the-wall radar 
imaging (TWRI) in both military and commercial 
applications, such as homeland security, urban counter- 
terrorism, search and rescue missions [1 ]-[7]. 


Through-the-wall radar images the targets behind the 
wall by transmitting ultrawideband (UWB) electromagnetic 
(EM) waves and processing the reflected signal from the 
wall and the targets. In Synthetic Aperture Radar (SAR) and 
Ground Penetrating Radar (GPR), the targets are usually 


stationary. However, in through-the-wall imaging the targets 
may move when the radar is collecting the scattered signal. 
During a long data collection time, the target may have 
moved such that the beamformed solution is no longer valid 
as the target is not in its previous position anymore. This will 
cause smearing, displacement and false targets in the image. 
Fast data acquisition avoids image smearing and imaged 
target displacements and allows prompt actionable 
intelligence in TWRI. In Sections II and II, we first present 
a generalized compressive sensing approach for reduced data 
collection and fast data acquisition in TWRI, and an efficient 
beamforming algorithm for through-the-wall imaging with 
multiple-input multiple-output (MIMO) radar which is real- 
time in data acquisition. Then, we present the detailed 
formulation and some representative results of a generalized 
2-D diffraction tomographic algorithm for real-time TWRI. 
In the following, we present some backgrounds on the above 
three topics. 


It is well known that achieving high imaging resolution 
in both cross-range and down-range requires the synthesis of 
a long aperture and transmission of an UWB signal. These 
requirements result in a large amount of space-time/space- 
frequency data and long data acquisition time in TWRI. 
Reduction of the data volume and acceleration of the data 
acquisition is important in TWRI, as it accelerates 
processing, and subsequently, allows prompt actionable 
intelligence [8]-[10]. It also relaxes constraints on system 
aperture and bandwidth, and creates different design and 
deployment paradigms, which are more flexible than those 
underlying conventional SAR operations. The capability of 
compressive sensing (CS) to reconstruct a sparse signal from 
far fewer non-adaptive measurements provides a new 
perspective to address these challenges in TWRI. In CS, it 
has been shown that a signal/image, which is sparse or has a 
sparse representation in some basis, can be captured from a 
small number of random non-adaptive linear projections onto 
the measurement basis. The original signal can be 
reconstructed with incomplete measurement data through 
convex optimization that uses the sparsity as an important a 
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priori information [11]-[13]. There have been several 
approaches that apply CS to radar imaging [14]-[18]. In [13], 
[14], [17], CS is applied to synthetic aperture radar to 
achieve the high resolution imaging with far fewer 
measurement. Gurbuz et al. proposed a novel CS data 
acquisition and imaging approach for step frequency and 
impulse ground penetrating radars in [15, 16]. In [8], [18] 
similar data acquisition and target reconstruction strategies 
were applied for stepped frequency continuous wave radar 
through-wall imaging in the framework of CS. In [9], CS is 
applied to impulse through-the-wall radar imaging. Other 
recent works on application of CS to TWRI are those 
reported in [19]-[20]. 

Most of the existing literatures on TWRI with CS deal 
with only a single layer wall, in which the ray tracing 
technique is used to calculate the wave propagation time as it 
travels through the wall. In many urban sensing applications, 
we often encounter situations to detect and identify targets 
inside a building with multilayered inner walls or walls 
separated by a hallway [3]-[5], [7]. This is challenging and is 
beyond the capability of existing through-the-wall 
beamformers using ray tracing techniques. In Section II of 
this paper, we expand on application of the CS to the 
problem of imaging of targets located behind either a single 
or multilayered building wall through the incorporation of a 
generalized forward model in the imaging algorithm to take 
into account and compensate for the wall effects. 


Existing TWRI are mainly discussed within the 
framework of monostatic SAR. One of the main drawbacks 
of monostatic radar system is the long data acquisition time 
to synthesize the aperture. It takes about a few minutes and 
over an hour to synthesize 1-D linear and 2-D rectangular 
apertures in the laboratory measurement, respectively [7], 
[22]. In Section HI, we address the problem of through-the- 
wall imaging with multiple-input multiple-output (MIMO) 
radar. The MIMO concept using multiple antennas at both 
the transmitter and receiver was first proposed to improve 
communication performance and later extended to radar 
systems and has been attracted increasing interest in TWRI 
in the recent years [25]-[29]. The MIMO radar system does 
not need to synthesize an aperture thus is real-time in data 
acquisition. Moreover, the spatially distributed multiple 
transmitters and receivers in MIMO radar increase the 
diversity of target illumination and scattering by viewing the 
target from multiple aspects [30]. 


In addition to fast data acquisition, the fast processing of 
the received signal to form an image of the targets behind the 
walls is critical to real-time monitoring of targets in TWRI. 
Up to present, several effective TWRI algorithms that take 
into account the wave reflection, bending, and delay effects 
due to the presence of the wall have been proposed, such as 
delay-and-sum (DS) beamforming algorithm, linear inverse 
scattering algorithms and time reversal based imaging 
algorithms, etc [6], [31]-[33]. Although successful imaging 
results can be obtained by using these algorithms, they are 
based on the pixel-by-pixel reconstruction of the image. The 
computation time increases significantly with the increasing 
of the number of pixels of the image, making these 
algorithms still not applicable to real-time processing, 


particularly when one deals with large scale imaging 
scenarios or 3-D imaging. TWRI algorithms must be 
computationally efficient, so that the imaging of the targets 
can be achieved within a few seconds/minutes for 2-D/3-D 
imaging with a portable computer. In [2], [34] 2-D and 3-D 
diffraction tomographic (DT) algorithms are presented for 
TWRI with a single layer homogeneous wall. The DT 
algorithm requires much less computational resources and is 
much more computationally efficient due to the easy 
implementation of the algorithm with Fast Fourier Transform 
and inverse Fast Fourier Transform (FFT/IFFT). In Section 
IV of this paper, we present a generalized diffraction 
tomographic algorithm for 2-D real-time TWRI that is 
capable of imaging of targets behind single or multilayered 
building walls. 


II. COMPRESSIVE SENSING FOR TWRI 
A. Forward model for imaging through layered walls 


Figure 1 shows a typical scenario of through-the-wall 
with monostatic SAR. The radar antenna moves along a scan 
line parallel to the homogeneous wall from Lmin to Lmax. The 
location of the m-th transceiver is denoted as rm. The radar 
operating frequency ranges from f min to f max with frequency 
step Af. The wall region may consist of single or 
multilayered walls. 


Under the point target model, the received scattered field 
can be written as 


E (tink, ) = I. dro(r)G(r,,.0.4,)G(r.4,,>, ) 


(1) 
where E, (Tp , k,,) is the received scattered field at the m-th 


antenna location, G (Im, r, kn) is the background layered 
media Green’s function, ø ( r ) is the reflectivity of the 
target, kn is the wave number of the n-th operating frequency. 





Lmin 


Fig.1 Measurement Configuration of TWRI 
Given the received signal at all antenna locations, the 
imaging function is given by [35] 
Lx CK 
max max —): 
I(r) =| J drp dk,E, (rnk, ) C7 (Cmr, k,) © 


where Amin and kmax are the wave-numbers of the minimum 
and maximum working frequency. It is noticed that in the 


p 
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above equation the reciprocity principle G (rm, r, kn) = G (r, 
Fim, An) has been used. To compute the image in (2), an 
efficient evaluation of the layered media Green’s function is 
required. This is complicated and time consuming as it 
generally involves the numerical or semi-analytical 
evaluation of the Sommerfeld integrals [36]. In order to 
compromise between the computation time and the exact 
evaluation of the Green’s function, we assume that the target 
is located in the far field of the antenna and approximate the 
far field Green’s function as [6], [7] 


ATR., 


where Rm is the distances from the m-th antenna location to 
the target and 7 is the transmission coefficients from the 
radar antenna to the target and will be given in Section II B. 





G(r 


tm ? 


r,k,)=T (3) 


Given the reflectivity of the target, by substituting (3) 
into (2) and ignoring the real-valued coefficient (47Rm) *, the 
received signal at the m-th receiver can be written as 


E (rnk )=| dro (rje? T? (r,r k, ) Ry ® 


tm? “n 2 tm? “n 


and the image can be then reconstructed as 


I(r)= [e úi dr „dk, E, (tk, ern - 





mın (5) 
=) 2 
T (r, En k, ) Rin 
B. Transmission through Multilayered Walls 
z 
Region 1 x E> By 
-d, 
Region 2 &, $ 
=d, 
SSS EEE 
Region N-1 
g a 
Region N Ex» Hn 


Fig.2 Transmission through Multilayered Media 


By incorporating the far field approximation of the 
layered media Green’s function, the imaging algorithm can 
be easily generalized to the imaging of targets located behind 
multilayered homogeneous walls, such as a brick wall, 
poured concrete wall or a hollow concrete wall. Let us 
consider the multilayered homogeneous media shown in 
Figure 2, where the permeability, permittivity and thickness 
of the n-th region are denoted as un, €n, and dn, respectively. 
Applying the boundary condition at the interface of each 
layer, the generalized reflection coefficient can be derived, as 
detailed in [23] 

+R 


i+1,i+2 


phe _ Rian t Rania @XP(J iss, (din —4)) (6) 
i,i+l 14R È -a exp(j2k,,.. (d, -d,)) 


i i+l i+ i 


where kiz is the normal component of wavenumber k; in the i- 
th layer, R., is the generalized reflection coefficient from 


the i+1 th layer to the i-th layer, and Ri;+; is the reflection 
coefficient from the i+1 th layer to the 7 th layer. The local 
half-space reflection and transmission coefficients are given 
as 


k_—k, | 
i = IZ JZ l T, = 2k, (7) 
Sk, +k, "ok, +k, 


The generalized transmission coefficient from the first 
layer to the N™ layer can also be derived from the continuous 
boundary condition at the interface of each layer [23] 


liN = exp(—jky.dy_, T jk,,d,)- 


(8) 
TT exp( jk,.(d, -dp )) S, pa 
=| 
where . 
T 
p-l,p rit (9) 


l= Kari exp( j2k,, (d, =d )) 
For a single layer homogeneous wall, the total reflection 
and transmission coefficients in (8) and (9) reduces to [6] 


R(1- 2k. (d,-—d 
g - Raltcexp (ih. (4-4) o 


I-(R,, j exp( 2k, (d, - d,)) 


(1-(R,) Jexp(j2h. (d,—4d,)) 
i a ae a a 
1-(R2) exp( j2k,, (d, —4,)) 


C. Imaging through multilayered building walls with CS 


Referring to Figure 1, we assume there are M receiving 
antenna locations and a transmitted stepped-frequency signal 
of N narrowband waveforms. The received signal at the m-th 
antenna with the n-th frequency bin is given by (4). Through 
discretization, (4) can be written in the following matrix 
form 


(12) 
where Zm is the measured data at the m ® antenna location 
with the 7 ® element [za] SE, (r 


tm? 


Z, = ¥,,8 


k,). s is a weighted 
indicator vector defining the target space, which is a K x L 
pixel image and vectorized into an K L x 1 column vector. 
The target space and the measured data is related to Py, 
whose (n, p) ® element is given by 
[P], p =T (estima k) eR 


p? tm? i (13 ) 

In contrast to conventional back-projection beamforming, 
which uses data measured at all antenna locations and all 
frequencies, in the CS through-the-wall imaging approach, 
we utilize a set of random measurements at a random subset 
of J frequencies at each antenna location, JU N . Then the 


measurement data at the m-th antenna can be written as 














(14) 
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where ®,, is a J x N measurement matrix in which each row 
has only one nonzero element, which is equal to one. The 
location of this nonzero element corresponds to the index of 
the measured frequency bin in the transmitting frequency 
sequence. In essence, the above measurement matrix is 
formed by randomly choosing J rows from Yn. 


For coherent processing of the received data, we 
superpose the matrix equation in (14) for all M illumination 
locations to form a large matrix 


y=OYV¥s=O0s (15) 


where 
y =[ylyloyt, | ’ P = diag|®,,®,,---®,, |. 


W=[PPOW wT) 00 


If the target space is sparse, which is generally valid for 
most TWRI applications, the reconstruction of s in the above 
equation can be achieved by solving the following sparse 
constraint optimization problem 


(16) 
For radar imaging problem, the measurement data are 


inevitably contaminated by noise. The data at the m-th 
antenna location can then be written in the following form, 


(17) 
where um is the measurement noise at the m-th antenna 
location. For noisy measurements, (15) can be rewritten as 


(18) 


S = argmin ||s||, s.t. y=Os 


Yn —~9,S8+u,, 


y=Ost+u 
p wie aT rY 
W ere u = u; „u, TEEPE. y . 


By separating the real and imaginary part of both sides of 
(15), it can be written in the following equivalent form 
y=I's (19) 


where 


y =| Re(y") Im(y")] . S= | Re(s’) Im(s’ ii 
Re(®) -Im(0) 
Im(®) Re(@) 
From [16], [24], robust reconstruction of a sparse signal/ 
image under noise corrupted data can be achieved by solving 


the following convex optimization problem which is also 
referred to as Dantzig selector 


S = argmin ||| s.t. | I’ (y -T8)| <ô 


T 


(20) 


where ò represents a small tolerance error which can be 
determined using the cross validation strategy for an 
automatic selection of the error value in the optimization 
process [21]. 


D. Numerical Results 


In the first example, we investigate the imaging of a 
rectangular PEC object behind a single layer homogeneous 
wall. The target is 0.3 m x 0.2 m and is 1.5 m away from the 
front boundary of the wall. The permittivity, conductibility 


and thickness of the wall are €&= 6, o = 0.05 S/m and d = 0.2 
m. The monostatic radar system scans the region of interest 
at a distance 0.3m away from the wall, synthesizing a 2 m 
aperture with 0.05 m inter element spacing. The operating 
frequency ranges from 1 GHz to 3 GHz with 61 equally 
spaced frequency bins. We have used background 
subtraction to remove the direct wall reflections [6]. 
Alternatively, in practice the effect of the wall reflections can 
be removed by either time-gating the received signal, or by 
analytically modeling the wall in frequency-domain and 
subtracting its response from the cumulative wall plus the 
target response [7]. 


Figure 3 (a) shows the imaging result using conventional 
back-projection beamformer in (5) using the full set of data, 
i.e. data measured at all frequencies and all antenna 
locations. For this single layer homogeneous wall, the 
transmission coefficient in (11) is used. If only two 
frequencies are randomly measured at each antenna location, 
imaging result using the conventional back-projection 
beamformer is shown in Fig. 3(a). Due to the lack of 
sufficient information about the target, the image is blurred, 
distorted, and has higher sidelobe levels, as evident from Fig. 
3(b). By exploiting the sparsity of the image space and 
solving the sparse constraint optimization problem in (20) 
with the same limited amount of measurement data, the 
reconstructed image is given in Fig. 3 (c). From the CS 
through-wall imaging result we find that the targets is clearly 
identified and well located at the true position. From the 
comparison of Fig. 3 (b) and (c) it is clear that a much higher 
resolution and cleaner image can be obtained through sparse 
constraint optimization. As seen from the CS imaging result 
with only 3.3% of the full set of data, the CS imaging 
approach gives even better and less cluttered image than that 
beamformed using full set of data in Fig. 3 (a). 
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(b) 


Down Range (m) 





0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 
Cross Range (m) 


(c) 
Fig. 3 Imaging of target behind a single layer wall (a) Beamforming result 
with full set of data; (b) Beamforming result with limited data; (c) CS 
imaging result with limited data. 


In urban sensing applications, we often encounter 
situations to detect and identify targets inside a building with 
multilayered inner walls or walls separated by a hallway. 
This is challenging and beyond the capability of the standard 
back-projection beamformers using ray tracing technique, 
and has not been studied in the existing CS TWRI 
approaches reported in the literature. In the second example, 
we present the imaging result of multiple targets behind 
external and interior walls separated by a hallway, as shown 
in Fig. 4 (a). The width of the hallway is 1.1 m, which is 
sandwiched by an exterior and internal wall whose dielectric 
constant and thickness are €, = 7 and d = 0.15 m. We note 
that the wall parameters in such scenario can be efficiently 
estimated using the technique described in [41-43]. Figure 4 
(b) shows the beamforming result of the two targets using the 
full set of data. If only two frequencies are measured at each 
antenna location, which results in 41x2 spatial-frequency 
measured data, it is very difficult to identify the targets as 
shown in Fig. 4 (c). However, a much less cluttered and 
focused image can be achieved by using the generalized CS 
TWRI approach with the same random measured limited 
data, as can be seen from Fig. 4 (d). 
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(d) 
Fig. 4 Imaging of target behind the inner wall (a) Simulation geometry; (b) 
Backprojection beamforming result with full set of data; (c) Backprojection 
beamforming result with limited data; (d) CS imaging result with limited 
data. 


IU. THROUGH-THE-WALL IMAGING WITH MIMO RADAR 


In this section, we present the detailed formulation of 2-D 
through-the-wall imaging with MIMO radar for real-time 
data acquisition. This method can be easily generalized to 3- 
D TWRI, if needed. 


A. Problem Formulation 


Figure 5 shows a 2-D scenario of through-the-wall 
imaging with MIMO radar. The MIMO radar system consists 
of M transmitting antennas located at [xm, Zm| and N 
receiving antennas located at [x Zm]. The operating 
frequency ranges from fmin tO fmax With a frequency step A f. 
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Targets 


Fig.5. Configuration of MIMO radar through-the-wall imaging 


Under the point target model, which ignores the multiple 
scattering effects, the received signal can be written as 


E, (Pins Pek) = [G (Pn P-k, )G(P: Pinsky )O(P) dP 


(21) 
where Es (Pin, Pim, Kp) is the received scattered field at the m- 
th receiver location due to the illumination of the n-th 
transmitter, o is the reflectivity of the target, pim, Pm and p 
are the position vectors of the m-th transmitter, n-th receiver 
and target, 1.€., Pim = (Xim, Zim), Prn = (Xrn, Zm), r = (x, Z), kp iS 
the freespace wavenumber of the p-th operating frequency, G 
(Pim, P, kp) and G (P, Pran, kp) are the layered media Green’ s 
functions, which relate the wave propagation process from 
the transmitter to the target and from the target to the 
receiver, in the presence of the wall. Then, the through-the- 
wall image can be reconstructed as 


I(p)= ji dk, f, “dp, f, ~ dP,E, (Pn Pmnskp) 


(G (Pn: Pk, )G(P Pinsky) 


In the above imaging formulation the inner two integrands 
corresponds to the beamforming over all the transmitter and 
receiver locations and the outer integrand is the coherent 
summation over all the operating frequencies. Assume that 
the targets are within the far field of the antenna, the far field 
Green’s function of the walls can be approximated as [6], [7] 


(22) 


— Tk y Rin 
G (P, Pinsky) =T, (0s Omsk, )—— 
p) 4r R, (23) 
Thy Run 
C (Pn: Pk.) =T. (Oo Proky GR 


where Rm = | Pm - P |, Rm = | Pm -p |, T; and T, are the wall 
transmission coefficients from the transmitter to target and 
the target to receiver, respectively. 


By substituting (23) into (22), the MMIO radar through- 
the-wall imaging formula can be derived as 





L 


I(p)= j dk, [ d Py - Ppt, (An Prnrk, ) 


T (P, Pask, ) T7 (Ps Pyk) E 


In the above equation, the real-valued term 167?RmRr is 
omitted. It is found from (24) that the wall effect is 
efficiently taken into account and well compensated through 
the background media Green’s function. Due to the 
incorporation of the layered media Green’s function, the 
imaging algorithm can easily be generalized for the imaging 
of targets behind multilayered building walls. 


(24) 


B. Numerical Results 


As the first example we investigate the imaging of targets 
behind a single layer homogeneous wall with MIMO radar. 
The forward scattering data was simulated using a full wave 
EM solver based on finite difference time domain method 
(FDTD). In the simulation we used eight transmitters equally 
spaced from -1 m to 1 m and thirty two receivers equally 
spaced from -1.25 m to 1.25 m along x direction at a standoff 
distance of 0.3 m. The dielectric constant, conductivity and 
thickness of the wall are £» = 6, op = 0.05 S/m, and d = 20 
cm. The operating frequency ranges from 1 GHz to 3 GHz. 
The three target investigated in this example are rectangular, 
cylindrical, and square PEC objects with the designation of 
Target 1, Target 2, and Target 3, respectively. The spatial 
characteristics of the three targets are summarized in Table 1. 


Table 1: Target Spatial characteristics 
Size 


Center Position 


Cylindrical R=0.1m (0, -1.9 m) 
Dmm | OSmeAFm 


In order to demonstrate the necessity and importance of 
taking into account the wall effect, first a conventional 
freespace beamformed imaging result is shown in Figure 6 
(a), where the true regions of the targets are indicated with 
white dashed rectangle, circle and square. From the freespace 
beamformed imaging result we find that the targets are not 
only shifted in the downrange (about 0.3 m) but also 
widened in the cross range. This is due to the fact that the 
wave propagation effect is not considered and compensated 
in the conventional freespace radar imaging. Figure 6 (b) 
shows the imaging results of the targets using the proposed 
MIMO through-the-wall beamformer. From this image we 
find that the three targets are clearly identified and well 
located at their true locations. It is clear from this image that, 
due to the proper compensation of the wall effects, high 
quality focused images of the targets can be achieved using 
the proposed beamformer for through-the-wall imaging with 
MIMO radar. 


In the second example, we present the imaging of target 
behind external and interior walls separated by a hallway 
using MIMO radar (as shown in Figure 7 (a)), which is a 
representative scenario in urban sensing of building interior 
targets [4], [7], [24], [37]. The permittivity, conductivity and 
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thickness of the exterior and interior walls are £» = 6, o» = 
0.05S/m, and d = 20 cm. The width of the hallway between 
the exterior and interior walls is 1.2 m. The targets under 
investigation are square and circular cylinders with the 
designation of Target 1 and Target 2, respectively. The 
spatial characteristics of the two targets are summarized in 
Table 2. The radar operating conditions are the same as in 
the previous example. 


Table 2: Target Spatial characteristics 





Cylindrical 





E 
N 
x(m) 
(a) 
E 
N 





x(m) 


(b) 
Fig. 6. Imaging of targets behind a single layer wall with MIMO radar (a) 
freespace imaging result; (b) through-the-wall imaging result. 


Figure 7 (b) shows the imaging result of the targets using 
freespace beamforming algorithm. The true regions of the 
targets are indicated by the white dashed circle and square. 
From the freespace imaging result we find that the targets are 
not correctly localized at their true locations. For multiple 
building walls, the targets are observed to be shifted about 
0.65 m away from their true locations in the down range, 
which is about twice the target dimension. Figure 7 (c) 
shows the imaging results of the targets using the proposed 
MIMO through-the- wall beamformer. From the image we 
find that the targets are clearly identified and well located at 
their true locations. It is clear from this figure that the 
proposed MIMO beamformer is successful in the imaging of 
the target behind multilayered building walls without 
distortion or displacement, and with a significantly reduced 


R=0.15 m (0.25 m, -2.95 m) 


data acquisition time as compared to conventional 
monostatic SAR operation. 
oO Tx 
O Rx 
o O O O O O 
1.2 m 


Hallway 





Interior Wall 





z(m) 








x(m) 


(b) 





z(m) 








“4 -0.5 0 0.5 1 
x(m) 


(c) 
Fig. 7. Imaging of targets behind external and interior walls separated by a 
hallway with MIMO radar (a) simulation geometry; (b) freespace imaging 
result; (c) through-the-wall imaging result. 


IV. GENERALIZED DT ALGORITHM FOR REAL-TIME 
TWRI 


In [34] a 2-D diffraction tomographic algorithm for the 
imaging of targets behind a single layer wall, derived from 
the stationary phase method, has been presented. In this 
section, we present a generalized DT algorithm for 2-D real- 
time TWRI based on the exploding reflection model that is 
capable of imaging of targets behind either single or 
multilayered building walls. 
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A. Problem Formulation 


As shown in the TWRI configuration with SAR in Figure 
1, the radar system moves along a scan line parallel to the 
wall in the observation region denoted as Doss. The radar 
operating frequency ranges from f min to f max. The wall 
region consists of either single or multilayered building 
walls. The targets to be imaged are located in an inaccessible 
investigated region denoted as Dim behind the walls. The 
investigated domain is chosen to be a rectangular region rom 
—Xiny tO Xiny and Zmin tO Zmax. The relative complex permittivity 
Er (r) inside Din, can be defined as 


e, (r)= l E, (r) r € Dp 


, (25) 
E,, otherwise 


where r = Xx + Zz is the target position vector. 


Let us assume that the electric current is a 2-D point 
source, which is equivalent to a 3-D line source directed in 
the y direction. Then the scattered electric field at the 
receiver can be written as 


E (rk) =i | G(r„r,k) E (r,r k)g(r)dr 6) 


2” tm? 


where r „= XX, + aie is the antenna location, G is the 


Green’s function for the background media, LE, (r r 4k ) is 


2“ tm? 


the total electrical field inside the target, x(r) is the 


contrast function of the object y (r) DE (r) —l,andé., is 


the relative permittivity of the object. Since the total field 
inside the target also depends on the contrast of the target, 
(26) is a nonlinear equation that is very time consuming to 
solve. In order to linearize the equation, the first-order Born 
approximation is employed, in which the total field inside the 
target is approximated by the incident electric field 


E, (r,r k) = E” (r,r ok) = iknyG (r,r k) D 


2 tm? 2 tm? 


Then the received scattered field can be derived as 
E (1,,k) = imk? | G(r,,..k)G(r,,.1.k) 7(r) dr (28) 


where 77, is the wave impedance in the freespace. 


Using the reciprocity property of the Green’s function 
G(r rk) = G(r r k), (28) can be rewritten as 


tm? 2“ tm? 


E (Prk)=ink |, G(PxP-k)z(p)dp 09) 


The Green’s function for the planar layered walls can be 
expressed in the spectral form as detailed in [21] 


T(K.) (ihe (Xr a 


Į o0 
T A 


1z 
(30) 
where 7, the transmission coefficient for the walls, can be 
calculated from (8). 


For derivational convenience, (29) can be written in the 
following operator form 


E,(1,,-«) =(Lz)(r) (31) 


The operator L maps the contrast function from the object 
space to the data space, in particular, 
































L(Q) imk | G(r," k)(dr (32) 
Then the image can then be expressed by 
I(r) =L'(E,(r,,,%)) (33) 


where Z° is the adjoint operator that maps the scattered field 
from the data space to the object space, in particular [38], 
[40] 
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Nu 





k * 
. max 3 2 
-imj Uke | [G(r tsk) | Oar, 

(34) 
where the * denotes the complex conjugation. The complex 
conjugation in the adjoint operator is essentially the back- 
propagation of the wave to the location where it is excited. 
From the concept of time reversal imaging, the complex 
conjugation in the frequency domain is equivalent to the time 
reversal in the time domain. The back-propagation of the 
time reversed field will focus on the target. 


In this paper, to accelerate the imaging algorithm the 
assumption of exploding reflection model is employed [39]. 
In particular, the squared term of the Green’s function 


G(r r,2k). Then (34) 


tm ? 


r,k ) is approximated by G (r 


tm ? 


can be further simplified as 
LOE -Ifi dk f, G(r at.28)( 


Substitute the spectral form of the layered media Green’s 
function (30) into (11) one has 


P(E (to) 


No i= kdk| dk, T (Hs) olkx-ikiszn tik yee 
Am -Fmin ag k 


J Dix, E, (Tim pk Je" dx 

















)dr,, (35) 


(36) 





lz 


It is noticed that due to the use of the exploding model, 


kee k Tal, ee 


formula can be derived as 
I(r) =L'(E,(r,,,*)) 


— o kmax idk | dk E, (k, k) 4 &) ote -ikizzR +iky,Z 


= i - 
Agr Krin 5 


Then the imaging 


(37) 
where E, (k.k ) is the Fourier transform of E, (r k) 


tm ? 


with respect to X,,, in particular 
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E,(k,.k)= I, E,(r,.k)exp(—ik,xp)dx, (38) 


The above integral for E, (k.k ) can be efficiently 
evaluated using FFT. 


B. Numerical Results 


In the first example, we present the imaging of a human 
behind a hollow concrete wall. The measurement data was 
generated using XFDTD®, a commercial full-wave 
electromagnetic simulator based on FDTD from Remcom 
Inc. The dimension of the HiFi male human model is 0.57 m 
x 0.324 mx 1.88 m and made up of 2.9 mm cubical FDTD 
mesh cells, 23 different tissue types with realistic dielectric 
and conductivity parameters. The front and side view of the 
human is shown in Fig. 8. The radar measures at a distance 
of 0.3m from the front wall synthesizing a 2 m aperture with 
inter-element spacing 0.05 m. In the XFDTD simulation it 
took about 1.5 minutes for each transmitting location on a P4 
8-core 2.8 GHz personal computer with 16G RAM memory. 
The operating frequency covers the range from 1-3 GHz. The 
inner and outer wall has the same electrical parameters with 
dielectric constant £r = 6 and conductivity o = 0.01 S/m. The 
thicknesses of outer and inner concrete walls are 0.1 m and 
are the same as the air gap between the two walls. 


Fig.8 Simulation geometry of a human behind a hollow concrete wall 
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Fig. 9 Imaging of the human behind a hollow concrete wall 


The investigation domain is a2 m x 2 m and is divided 
into 100 x 100 pixels. Figure 9 is the imaging result of the 
human using the proposed imaging algorithm. The 
approximate true region of the human is indicated with a 
white dashed ellipsoid. From this figure we find that through 
the proper incorporation of the multilayered media Green’s 
function, a high quality focused imaging result can be 


achieved using our proposed imaging algorithm and is very 
computationally efficient. It took only about 0.63 s in Matlab 
to form the image of the human on a P4 4-core 2.3 GHz CPU 
personal computer with 4G RAM memory. It can be 50-100 
times faster than existing imaging algorithms in [6, 31 ,32]. 

In the last example we investigate the problem of imaging 
of multiple targets located behind external and internal walls 
separated by a hallway, which is a typical scenario in urban 
sensing of building interior targets. The thickness of both 
the external and internal walls is 0.2 m with relative 
permittivity € , = 6 and conductivity o= 0.01 S/m. The width 
of the hallway between the two walls is 1.2 m, as shown in 
Figure 10 (a). The radar antenna transmits an UWB EM 
wave covering the frequency range from 1GHz to 3GHz and 
measures the return signal at a distance of 0.3 m from the 
external wall over a synthetic aperture from -2 m to 2 m. The 
two targets investigated in this example are two PEC 
cylinders with a radius of 0.3 m and 0.2 m centred at (-0.7 m, 
-3.6 m) and (0.8 m, -3.3 m), respectively. The investigation 
domain is a 2.4 m x 2 m rectangular region and is divided 
into 120 x 100 pixels. 
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Fig.10 Imaging of multiple targets behind the external and interior walls 
separated by a hallway (a) simulation configuration; (b) imaging result 


Figures 10 (b) shows the imaging results of the targets 
using the proposed algorithm, where the correct locations 
of the targets are indicated with white dashed circles. From 
Fig. 10 (b) we find that the two targets are clearly identified 
and well localized at their correct locations. It is evident 
from the image that, the propagation effect has been 
properly taken into account and high quality focused 
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images of the targets can be achieved using the proposed 
imaging algorithm. In this example, it takes only 0.67 s to 
form the image shown in Fig. 10 (b) when the generalized 
DT TWRI algorithm is used. A considerable time saving 
has been achieved in the process of imaging due to the 
utilization of FFT in the algorithm. 


V. CONCLUSIONS 


Conventional through-the-wall radar imaging techniques 
face considerable challenges, such as requirement to acquire 
a large amount of measured data, long data acquisition and 
lengthy processing time. In this paper, we present three 
techniques to address these challenges in TWRI: 1) 
compressive sensing imaging approach; 2) through-the-wall 
imaging with MIMO radar; 3) generalized DT algorithm for 
real-time TWRI. All these techniques are applicable to the 
imaging of targets behind either a single or multilayered 
building wall via the incorporation of the layered medium 
Green’s function that fully take into account and compensate 
for the wall effects. Representative results show that the 
proposed techniques are very efficient to reduce the long 
data acquisition and processing time in TWRI and can 
provide high quality focused imaging results in various wall- 
target scenarios. All these techniques can be generalized to 
3-D TWRI scenarios and the results will appear in the future. 
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